Breakdown of the Fermi-liquid regime in the 2D Hubbard model from a two-loop 
field-theoretical renormalization group approach 
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We analyze the particle-hole symmetric two-dimensional Hubbard model on a square lattice start- 
ing from weak-to-moderate couplings by means of the field-theoretical renormalization group (RG) 
approach up to two-loop order. This method is essential in order to evaluate the effect of the 
momentum-resolved anomalous dimension r?(p) which arises in the normal phase of this model on 
the corresponding low-energy single-particle excitations. As a result, we find important indications 
pointing to the existence of a non-Fermi liquid (NFL) regime at temperature T — > displaying a 
truncated Fermi surface (FS) for a doping range exactly in between the well-known antiferromagnetic 
insulating and the d x 2 _ H 2-wave singlet superconducting phases. This NFL evolves as a function of 
doping into a correlated metal with a large FS before the d x 2_ y 2-wa,ve pairing susceptibility finally 
produces the dominant instability in the low-energy limit. 

PACS numbers: 71.10.Hf, 71.10.Pm, 71.27.+a 



I. INTRODUCTION 



The physical nature of underdoped cuprates both 
above and below the superconducting temperature T c 
is still subjected to strong debate 1 -^! 3 -. The recent ex- 
periment of Doiron-Leyraud et olA applying a mag- 
netic field strong enough to destroy the superconducting 
state established the existence of quantum oscillations 
and coherent Fermi surface (FS) pockets for this regime. 
By contrast, earlier angle-resolved photoemission exper- 
iments (ARPES) reported only observations of discon- 
nected Fermi arcs located in the nodal regions of mo- 
mentum spaced. Moreover, according to both scenarios, 
there exist charge pseudogaps and no quasiparticlc-likc 
excitations in the corresponding antinodal sectors cen- 
tered around (±7r,0) and (0,±7r). The absence of these 
quasiparticle excitations is signalled by the vanishing of 
the quasiparticle peak Z(p) in the spectral function at 
these antinodal regions. The pseudogap behavior also 
reflects itself in the underdoped superconducting phase. 
One manifestation of this is the fact that the supercon- 
ducting gap and the pseudogap behave distinctively as 
a function of doping: the superconducting gap located 
in the nodal regions centered around (±7r/2, ±7t/2) de- 
creases continuously whereas the pseudogap continues to 
increase as the hole doping is reduced from its optimal 
valued. Such a nodal-antinodal dichotomy is therefore 
a common trend in underdoped cuprates. 

A minimal model for describing the dynamics of 
the interacting electrons in these high-T c superconduc- 
tors is the single-band 2D Hubbard model 9 . This 
model has been investigated extensively with sev- 
eral numerical techniques such as exact diagonal- 
ization, quantum Monte Carlo and the more re- 
cent quantum clusters approache o 10 ' 11 , as well as by 



semi-analytical functional renormalization group (RG) 
methodsi 2 .^^^^ 1 ^ 1 ^^. Several of these works 
successfully obtain an antiferromagnetic phase in the 
model near half-filling and the onset of a d x 2_ y 2- 
wave singlet superconducting phase away from half- 
filling when the temperature is lowered below a critical 
valu o 10 i 12 ' 13 i 14 . These results therefore give further sup- 
port to the point of view that the 2D Hubbard model 
might indeed capture the essential aspects of the physics 
displayed by those strongly-correlated materials. 

On the experimental side, the insulating antiferromag- 
netic phase in the hole-doped cuprates is quickly de- 
stroyed at a very small but nevertheless nonzero dop- 
ing. If we approach this phase from a larger doping 
regime, this suggests that the existing Fermi surface 
should shrink to zero with the corresponding quasiparti- 
cle weight Z(p) becoming suppressed all along the under- 
lying FS at a possible quantum critical point. It is there- 
fore clearly important to account for the physical effects 
associated with such a vanishing of Z at a nonzero dop- 
ing near half-filling. This analysis has been initiated in 
recent years for the 2D Hubbard model and its extensions 
in the context of the fermionic RG framewor k 21 i 22 i 23 i 24 . 
However, since the self-energy effects manifest themselves 
only at two-loop order or beyond, in order to devise con- 
serving many-body approximations 2 ^ for this model, one 
should necessarily evaluate all flow equations up to the 
same order of perturbation theory. For this reason, we 
implement in this work a two-loop field-theoretical RG 
calculation as a function of doping concentration for some 
important quantities in the 2D Hubbard model. 

Following previous one-loop RG calculations per- 
formed in the 2D Hubbard model, we choose to consider 
here only weak-to-moderate initial couplings. As a result, 
we find important indications pointing to the existence 
of a non-Fermi liquid (NFL) regime at a nonzero doping, 
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FIG. 1: The g-ology parametrization for the marginally rele- 
vant interaction processes of the 2D Hubbard model, using as 
a reference the corresponding noninteracting FS at half-filling 
for simplicity. 



before the d x 2_ y 2--wave singlet superconducting instabil- 
ity finally dominates over the antiferromagnetic fluctua- 
tions. An important point we wish to stress here is that 
this NFL regime emerges from the nonzero momentum- 
resolved anomalous dimension rj(p) contribution to the 
single-particle excitations, which arises only at a two-loop 
RG level or beyond close to half-filling. In the corre- 
sponding one-loop RG flow equations, even if the initial 
values are taken to be reasonably small, the renormal- 
ized interactions rapidly flow towards a strong coupling 
regime in the low-energy limit. This fact also indicates 
the importance of higher-order quantum corrections to 
the full description of the low-energy dynamics of the 
2D Hubbard model. Our work therefore takes seriously 
this observation and, for this reason, it represents a step 
forward in this direction. 

This paper is organized as follows. In Sec. II, we 
explain the methodology employed to discuss the 2D 
Hubbard model. In this part, we will choose to be 
very schematic since the fermionic field-theoretical RG 
methodology was already explained at length and in full 
detail in the context of a simpler 2D flat FS model 
elsewhere^. Our main emphasis will be rather to high- 
light the final integro-differential two-loop RG flow equa- 
tions resulting from the application of this method to the 
2D Hubbard model. In Sec. Ill, we move on to the nu- 
merical solution of these RG flow equations. Lastly, in 
Sec. IV, we present our final conclusions regarding our 
two-loop RG calculation and we point out open issues 
that still have to addressed for the full clarification of 
this important problem. 



II. METHODOLOGY 

We start by defining the Hamiltonian of the 2D Hub- 
bard model in momentum space 



k.cr \ »/ Pj k,q 

(2.1) 

where ipL and t/w ar e the usual fermionic creation 
and annihilation operators with momentum k and spin 
projection a, t is the electronic hopping amplitude to 
nearest neighbor sites, fj, is the chemical potential which 
controls the band filling (and, consequently, the doping 
parameter), U is the local on-site repulsive interaction 
and N s is the total number of lattice sites. Since we 
will be interested only in the universal quantities of this 
model, we can linearize the tight-binding energy disper- 
sion £k = —2t\ L cos{k x ) + cos(fcj,)l — fi around the FS as 
£k ~ Uj?(k)n.(k— kp(/i)) with the Fermi velocity given 

by v F (k) = |(Vk£k|k=k F 0*))| = 2£yW* k x + sin 2 k y , 
kp (/i) being the Fermi momentum which defines the non- 
interacting FS for a continuous doping parameter and n 
is a unit vector perpendicular to the FS. 

The thermodynamical properties of this model can be 
computed from the coherent-state Grassmann represen- 
tation of the partition function 

Z = J V$,^]e-Io dT ( Lo ®^ +Li ^'^\ (2.2) 

where /3 — 1/T. The noninteracting Lagrangian is de- 
fined in a standard way 

LohM = Z / MKr){d T +^)MKr), (2.3) 

where J" k = J ra^p and the interacting Lagrangian in 
turn reads 

L mt $,ip} / / g(pi,p 2 ,P3) 

a a , JPI JP2 JP3 

XiMPi +P2 - P3,t)^ ct /(p 3 ,t)?/v(p2,t)Vv(pi,t). 

(2.4) 

The Eqs. (|2.3p and (|2.4p therefore define our bare quan- 
tum field theory which is regularized in the ultraviolet by 
restricting the momenta to |k| < Ao, where the cutoff is 
chosen to be Ao = At. Since this should correspond orig- 
inally to the 2D Hubbard model, we must set the bare 
interaction g{pi, P2, P3) initially equal to the local inter- 
action U. However, as we will see next, this functional 
dependence of the coupling on the momenta will in fact 
play a crucial role in the low-energy effective theory. 

Here we perform all the calculations in the T — > 
limit. The starting point of our approach is the nonin- 
teracting FS of the 2D Hubbard model defined by kp (n) 
for a continuous doping parameter. To include the effect 
of interactions, we use a g-ology parametrization adapted 
appropriately to our 2D problem (for more details on this 



procedure in the context of a simpler 2D flat FS model, 
see, e.g., Ref<2£). By means of a power counting analysis, 
one can easily verify that the dependence of the coupling 
function g (pi, p 2 , P3) on the components of the momenta 
normal to the FS are irrelevant in the RG sense. There- 
fore, we are allowed to project these coupling functions 
on the FS in such a way that their only functional de- 
pendence will come from the components parallel to the 
FS of the three external momenta pi, p 2 , and p 3 (the 
fourth momentum is given naturally by momentum con- 
servation p 4 = p 1 + p 2 — P3). In this way, we consider 
here the interaction processes which lead to singularities 
within perturbation theory in the low-energy limit. They 
are shown schematically in Fig. 1. These processes are 
the marginally relevant couplings in our RG theory, since 
their contributions become increasingly important at low 
energies. We neglect here the so-called (74-processes of in- 
teracting particles belonging to the same FS sector. Ex- 
perience with one-dimensional systems suggests that they 
should not alter qualitatively our results. Besides, we also 
neglect all marginally irrelevant interactions, since these 
contributions are not expected to change the universal 
properties of the model. 
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FIG. 2: The self-energy diagrams up to two-loop order. The 
one-loop term is the Hartree diagram and the two-loop con- 
tribution is the so-called sunset diagram. 
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(2.6) 



The methodology of our RG scheme follows closely the 
field-theoretical metho d 2 L 2 £ . In perturbation theory, in- 
frared logarithmic singularities typically emerge in the 
low-energy limit at the calculation of several quantities 
in the model such as the vertex corrections, the quasi- 
particle weight and the susceptibilities. We circumvent 
this problem by introducing appropriate counterterms 
at a flowing RG scale parameter A, in such a way that 
all the observables - i.e. the renormalized quantities of 
the theory - remain well-defined in the low-energy limit 
(A — * 0). As we have already pointed out in an earlier 
paper 26 , the problem of formulating a RG theory asso- 
ciated with a FS of a 2D model often requires the defi- 
nition of counterterms which are continuous functions of 
the parallel momenta along the FS. In this way, we must 
perform at two-loop level the following substitutions for 
the fcrmionic fields 



Vv(p,t) -» ^ /2 (p)Vv(p,T) 



Vv(p,t) 



Z A /2 (p)^(p,r), 



(2.5) 



where, from now on, all the momenta will correspond to 
the momentum projected on the FS. Besides, Z\(p) is 
the RG flowing momentum-resolved quasiparticle weight 
and it is naturally related in the limit of A — ► to the con- 
ventional many-body definition of the quasiparticle peak 
Z{p) = (l-dReE(uj,p)/doj\ u)=0 )- 1 . Using Eq. (|2"3]l . we 
also find that the bare and renormalized coupling func- 
tions are related to each order by 



where i =1, 2, 3, 3X, and BCS. The renormalized quanti- 
ties (labeled by the subscript R) generally depend on the 
RG scale A. In contrast, all the bare quantities will be de- 
noted here without any additional index. The functions 
&g] l R OP (pi , P2 , Pa ; A) and Agf™ ps (pi , p 2 , p 3 ; A) repre- 
sent the counterterms necessary to regularize at one-loop 
and two-loops respectively the four-point one-particle ir- 
reducible (1PI) vertex corrections in each of the corre- 
sponding parametrized scattering channels. 

Following the same strategy as it was described in full 
detail in Ref. 26 , we now calculate the quasiparticle weight 
Z\(p) at two-loop level starting from the definition of the 
renormalized self-energy of the present model. The cor- 
responding Feynman diagrams up to two-loop order are 
shown in Fig. 2. The one-loop diagram (i.e., the Hartree 
term) is generally independent of the external frequency 
lo. As a result, this contribution never renormalizes the 
quasiparticle weight. In fact, it only generates a constant 
shift in the chemical potential of the model which must 
be appropriately subtracted by a counterterm in such a 
way that the density of particles in the system always 
remains fixed during the RG flow. By contrast, the two- 
loop contribution (i.e., the sunset diagram) is the first 
contribution to the self-energy which produces a nonan- 
alyticity as a function of the external frequency lo. For 
this reason, this term alone will be responsible for the 
renormalization of the quasiparticle weight Za (p) at this 
level of perturbation theory as we vary the RG scale A 
towards the low-energy limit. 

The momentum-resolved anomalous dimension is con- 
ventionally defined by rj(p) = Ad In Za(p) /dA. From this 
expression, we get 



r?(p) = i 



dk / dq 

FS JFS 



A 



A + 2|/i|F(p,k,q) 
x [2 5 i fl (p,k + q-p,k) ff i fl (k,q,p) 
+ 2.g 2 fl(p,k+q-p,q)52ii(k,q,k + q-p) 

- gifl(P,k + q-p,k) ff2 fl(k,q,k + q-p) 

- g 2 fl(p, k + q p, q)sifl(k, q, p) + 2g 3R (k, p, q) 
x g 3 fl(k, p, q) - g 3R (k, p, q)s3ii(p, k, q)] 

1 \ / 1 



v F (k) +v F (q) J \v F (q) + v F (k + q-p) 



(2.7) 

where the integrals over the momenta now are sim- 
ply along the curve defined by the FS and, besides, 
F(p, k, q) = vp(k + q — p)/vp(k). It is interesting to 
note here that if we take the ID limit of the above equa- 
tion, we reproduce exactly the well-known result for the 
anomalous dimension of the Luttinger liquid at two-loop 
order as was first calculated long ago in Rei~. 

To derive the two-loop RG flow equations for the renor- 
malized coupling functions, we must take into account 
the fact that the bare theory does not know anything 
about the RG scale A. In other words, the bare couplings 
do not depend on this scale, i.e. Ad<?i(pi, P2, pa)/dA = 
0. Therefore, using Eq. (|2.6[) . we finally obtain 



d 1 4 

A^-5 , J i?(pi,P2,P3) = - 53 r ?(P I )fifl(Pl>P2,P3) 

t— 1 

-^j^^9m° P (Pi, P2, P 3 ) - A— A^^p^pa, p 3 ), 

(2.8) 

where i =1, 2, 3, 3X, and BCS. The Feynman diagrams 
for the vertex corrections necessary for the determination 
of the counterterm functions arc shown schematically 
in Fig. 3. These RG flow equations for the renormal- 
ized couplings are in fact complicated integro-diffcrcntial 
equations coupled to one another. We will discuss their 
numerical solution in the next section. 

To investigate what are the enhanced correlations in 
low-energy limit of the 2D Hubbard model, we must cal- 
culate the linear response of the system to small external 
fields. Therefore, we must add to our original Lagrangian 
the folllowing term 




FIG. 3: Some Feynman diagrams showing the vertex correc- 
tions up to two-loop order. 



superconducting and density-wave orders respectively^ . 
This added term will generate new Feynman diagrams - 
the three-legged vertices displayed in Fig. 4 - which will 
also generate new logarithmic singularities in the low- 
energy limit of our quantum field theory. Therefore, we 
must regularize these divergences by defining new coun- 
terterms as follows 



7 ^(p,q) = V^(p)VVp + q) 



-1/2, 



-1/2, 



7flfsc(P'<i; A ) 



a/3 / 



Z^ 1/2 (p + q )Z^ 1/2 (p) 



+ A 7 ^(p,q;A) 



7fl5w(p><i; A ) 



(2.10) 



where the Z-factors, as before, come from the redefinition 
of the fermionic fields at two-loop RG level displayed in 
Eq. (|2.5p . As a result, using again the fact that the bare 
response vertices do not depend on the RG scale A, we 
finally get 



! [^o(q)7so(P.«0^«(P.r)^(-p + q,r) A ^tow(P^) = 5 [f(P + <l) + V(p)}j%%w(p> <1 



-[r](p) +??(-p + q)]7flfsc(P J q) 



A^A 7 ^ sc (p,q) 



+ ''wfqlTwtp.qlifp + q^l^tp.T) +h.c. 

(2.9) 

where hsc(<l) and frow(q) a re the external fields and 
7sc(P'l) an d 7zuV(P' < l) are tn e response vertices for 



(2.11) 



Now, if we symmetrize the response vertices 75^(p,q) 

and 7j3v^(P)q) with respect to their spin indices we ob- 
tain the following order parameters 



5 



7ssc(p, q) = 7sc(P> q) - 7sc(P> q) 

7tsc(p, q) = 7sc(p j q) + 7sc(P. q) 
7sdw(p, q) = Jd W (p, q) - 7zw(p> q) 
icdw(p, q) = 1d W (p, q) + 7z»V(p> 



(2.12) 



where SSC and TSC correspond to singlet and triplet 
superconductivity and SDW and CDW stand for charge 
and spin density waves, respectively. Hence, using the 
above relations one can readily derive the RG flow equa- 
tions for each response vertex associated with a potential 
instability of the normal state towards a given ordered 
(i.e. symmetry-broken) phase. The initial conditions for 
these RG flow equations will determine the symmetry of 
the order parameter. Thus 



Xssc/rsc Cf) 




FIG. 5: Feynman diagrams for the susceptibilities associated 
with both singlet and triplet superconductivity and charge 
and spin density wave orders. 



charge and spin density wave susceptibilities. In addi- 
tion, the vector Q = (n, n) is the commensurate antifer- 
romagnetic wave vector. 



Hr(p, q; A = A ) = 1 (s-wave) 

7ir(p, q;A = A ) = -^=(cosp x - cosp y ) (d x 2_ y 2-w&ve) 



V2" 



(2.13) 



where i = SSC, TSC, CDW and SDW. Once we com- 
puted the response vertices associated with these order 
parameters, we can now proceed to calculate their cor- 
responding static susceptibilities. From Fig. 5, one can 
easily verify that they are given by 

A-^-Xfl,sc(q = 0;A) = - f -^[7i?,sc(p,~P;A)] 2 
dA J v F (p) 

A^x,.w(q = Q; A) = - (^^) / 

xfeflfffp + Q,p; A)] 2 

(2.14) 

where, to avoid cluttering up the notation too much, we 
simply write SC for both singlet and triplet and DW for 
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FIG. 4: The Feynman diagrams for the three-legged response 
vertices associated with superconducting and density-wave or- 
ders. 



III. NUMERICAL RESULTS 

To solve all these integro-differential RG equations nu- 
merically, we will follow the standard procedure of dis- 
cretizing the full FS continuum into N patches and then 
apply the fourth-order Runge-Kutta method. All numer- 
ical results in this paper will be presented for N = 32. 
Despite this moderate choice of number of patches, we 
point out that our results show good convergence prop- 
erties in the low-energy limit. 

First, we focus our attention on the numerical solu- 
tion of Eq. (|2.8[) for the renormalized couplings as a 
function of the ratio A/4t. As an initial condition for 
these equations, we choose the weak-to-moderate bare 
interaction U — 3t. Our results are displayed in Fig. 6 
for the case of /i = — 0.05£, which corresponds to a hole 
doping x of approximately x w 3.4%. In this regime, 
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FIG. 6: (Color online) The two-loop RG flow of the renor- 
malized couplings (in units of t) for initial interaction U = 3t 
and \i = — 0.05£. We discretize the FS into N — 32 patches, 
where 0, 8, 16, 24 represent the antinodal directions, and 4, 
12, 20, 28 refer to the nodal points. 
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FIG. 7: (Color online) The renormalized quasiparticle weight 
ZaQl) in the nodal and antinodal directions versus A/4/ for 
initial interaction U = 3t and initial condition ZA =4t(k) = 1: 
a) \n\/t = 0.045, b) \y\/t = 0.05, c) \fi\/t = 0.08 and d) 
\fj,\/t = 0.1. The solid lines represent curve fittings of the RG 
data. 



we observe that the inclusion of two-loop quantum fluc- 
tuations has important consequences to the flow of all 
couplings. Instead of exhibiting a divergence at a finite 
energy scale as obtained in several earlier one-loop RG 
investigation s 12 ' 13 ! 14 , the couplings now display a ten- 
dency to level off in the low-energy limit. This result 
therefore implies that the two-loop RG corrections clearly 
diminish the importance of fluctuations as compared to 
the one-loop RG theory. It is true, however, that several 
(but not all) couplings become saturated at fairly strong- 
coupling hxed values. This is a well-known problem in 
RG theory and happens as well in other applications of 
this method to quantum field theories in which fluctua- 
tion effects arc known to be very strong (the most no- 
torious example being the Wilson-Fisher fixed poin t 31 ' 32 
in </> 4 -theory at three dimensions). Nevertheless, the RG 
results in most cases are qualitatively correct even in a 
strong-coupling regime. The same is true in our case. 
When one of the renormalization couplings reaches an 
upper bound at the FS the whole flow is automatically 
stopped in spite of the fact that all couplings will not 
behave in the same way. We test the validity of our 
results by analyzing several physical quantities simulta- 
neously. Until this critical RG scale is reached by one of 
the couplings, no unphysical trend is observed in our cal- 
culations. In view of this, we hope that our two- loop RG 
results will also capture at least qualitatively the most 
essential aspects of the physics of the 2D Hubbard model 
in the low-energy limit. 

Next we present the results for the two-loop RG flow of 
the associated momentum-resolved quasiparticle weight 
Z\(p). Here Z\(p) is determined for several doping 
regimes as the low-energy limit is approached. The RG 
flows are displayed in Fig. 7 for the special choices of 
momentum at the nodal [near (7r/2,7r/2)] and antinodal 
[near (7r, 0)] directions. For doping levels up to a criti- 



cal value \i c = —0.045/, corresponding to x c « 3%, we 
observe that the quasiparticle weights associated with 
both nodal and antinodal directions become strongly sup- 
pressed and eventually scale down to zero at sufficiently 
low-energy scales. This result in fact holds for all points 
of the underlying FS of the system. From the two-loop 
RG flows of the several order-parameter susceptibilities 
displayed in Figs. 8(a) and 8(b), we observe that the 
s-SDW is by far the dominant instability at these low 
doping values. From these two results we can infer that, 
for \fi\/t < 0.045, the Fermi surface (if we define it sim- 
ply as the locus in k-space which exhibits a discontin- 
uous jump in the momentum distribution function) is 
completely smeared out by interactions and the resulting 
ground state should be given by an antiferromagnetic 
insulator. This agrees for instance with the fact that, 
exactly at half-filling, this model is known to have an 
antiferromagnetic insulating ground state for any finite 
local interaction U. Our two-loop RG approach there- 
fore successfully reproduces this important aspect of the 
problem. 

Following this, we probe the system at a slightly larger 
hole doping regime, i.e. fi/t = —0.05, which corresponds 
to x ~ 3.4%. At this regime, we obtain that the quasipar- 
ticle gap immediately closes around the nodal directions 
but remains unaffected in the antinodal regions. The 
complete suppression of the spectral weight only in the 
antinodal zones suggests the existence of a NFL phase 
with quasiparticle-like excitations only around the nodal 
directions of the Fermi surface. From a technical point 
of view, this result originates from the anisotropy in the 
anomalous dimension ?7(p) displayed in Eq. (|2.7p . This 
effect is to some extent reminiscent of the nodal-antinodal 
dichotomy and is in agreement with the general observa- 
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FIG. 8: (Color online) The two-loop RG flows of several 
static order-parameter susceptibilities for initial interaction 
(7 = 3/ with: a) and b) \n\/t = 0.045, c) \fi\/t = 0.08 and 
d) \/J,\/t = 0.1; CDW and SDW stand for charge and spin 
density waves; SSC and TSC stand for singlet and triplet su- 
perconductivity; the prefixes "s" and "d" refer to s-wave and 
d x 2_ y 2-wsve symmetries. 
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tion that, at light dopings, there are no coherent quasi- 
particle peaks at the antinodal points in the pseudogap 
phase of underdoped cuprates. However, we point out 
that the anisotropy obtained here in this work is still 
very small as compared with the experimental values ob- 
served in these materials. This is undoubtedly related to 
our choice of initial interaction U = 3t. If we increase 
the on-site interaction of the model, we expect that this 
trend will be amplified further. This conjecture finds sup- 
port in recent numerical results obtained from quantum 
cluster approaches 1 - 1 - for U — 8t, where a stronger nodal- 
antinodal dichotomy is indeed observed in the model at 
finite doping. 

Such a NFL character is observed here until the cor- 
responding doping reaches \n\/t = 0.08 (i.e. x » 5.1%). 
At this hole doping, both the antinodal and nodal quasi- 
particle weights become nonzero even when we extrap- 
olate our RG flow to very low A/4t scales. This indi- 
cates that the resulting FS is fully reconstructed in this 
regime. Moreover, as shown in Fig. 8(c), we observe 
already the beginning of the unlimited growth of the 
d x 2 _ y2 -wave pairing instability in the presence of the still 
leading s-SDW susceptibility. This mixed d-SSC and s- 
SDW metallic state is observed until \fi\/t = 0.1 (i.e. 
x sa 6.2%). At this doping value, as indicated in Fig. 
8(d), the d-SSC susceptibility finally overcomes the s- 
SDW and the system turns into a d rE 2_ J) 2-wave supercon- 
ductor. As a consequence, we point out that this latter 
result adds further support to the possible existence of 
a d x i_ y 2-w&ve superconducting state in the 2D Hubbard 
model as was first obtained in Refsj 12 ' 13 i 14 within a one- 
loop RG scheme. 

One important point we would like to emphasize here is 
that, for doping regimes satisfying < 0.08f , the domi- 
nant fluctuations in our results are always of s-SDW type 
(i.e. antifcrromagnetism) , in qualitative agreement with 
several previous studies focusing on this mode l 10 ' 12 i 13 i 14 . 
This result therefore favors the interpretation that the 
driving mechanism underlying our NFL evidence should 
be mediated by antiferromagnetic spin- fluctuations. 

We summarize our results displaying in Fig. 9 a two- 
loop RG phase diagram of the 2D Hubbard model as a 
function of both a critical RG scale A c and the hole dop- 
ing represented by fx. This phase diagram however must 
be interpreted only qualitatively since the numerical de- 
termination of the exact critical scale where the order- 
parameter susceptibilities truly diverge is of course very 
difficult and, for this reason, this is always based on a 
chosen criterion. Therefore, as we saw before, we find ev- 
idence of an extended antiferromagnetic insulating (AFI) 
phase from half-filling up to \n\/t = 0.045. For 0.045 < 
\fi\/t < 0.08, we obtain instead that the formation of a 
NFL metallic phase with a truncated Fermi surface in the 
nodal regions is favored. For 0.08 < \fi\/t < 0.1, the FS 
is fully restored with the resulting metallic state display- 
ing enhanced antiferromagnetic and d x 2 _ y2 -wave pairing 
correlations. Finally, for \fx\/t > 0.1, the system shows 
a d x 2 _ y 2-~w&ve superconducting phase. Moreover, if we 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 



FIG. 9: (Color online) Two-loop RG phase diagram for the 2D 
Hubbard model as a function of the critical RG scale A c and 
the hole doping \i for local on-site bare interaction U = 3t. 

suppose that it is possible to associate the RG scale A c 
with an effective critical temperature scale T c , this could 
in principle provide approximate estimates of the critical 
temperatures for the various phase transitions displayed 
by this model obtained within our two-loop RG scheme. 
Following this approach, we obtain that, for the d x 2_ y2 - 
wave superconducting phase, the highest effective critical 
temperature scale displayed in our two-loop RG phase di- 
agram is given by the curve fitting value A c /4i sa 0.004< 
(see inset), which is surprisingly in qualitative agreement 
with the result T c w 0.023i obtained in the literature for 
U = At within quantum cluster methods^. 



IV. CONCLUSION 

We have analyzed the 2D Hubbard model starting from 
weak-to-moderate couplings by implementing a func- 
tional generalization of the field-theoretical renormaliza- 
tion group approach up to two- loop order. Our calcula- 
tions here were mainly restricted to the limit of T — > 0, 
i.e. to the analysis of the universal ground state proper- 
ties and its low-lying excitations, as a function of doping. 
The two-loop RG scheme was necessary in order to dis- 
cuss the effect of the momentum-resolved anomalous di- 
mension 7y(p), which shows up in the normal phase of the 
model on the corresponding low-energy single-particle ex- 
citations. As a result, we have found evidence for an ex- 
tended antiferromagnetic insulating (AFI) ground state 
in the vicinity of half-filling. This result was particularly 
encouraging since it agrees with the fact that, exactly 
at half-filling, this model is known to produce an AFI 
ground state for any finite local interaction U. Our two- 
loop RG approach therefore successfully reproduced this 
important aspect of the problem. 

For a slightly higher doping regime, we have obtained 
that the formation of a NFL metallic phase with a 
truncated Fermi surface in the nodal regions was fa- 
vored. This effect is originated by the anisotropy of the 
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momentum-resolved anomalous dimension ?y(p) at this 
regime and it is reminiscent, to some extent, of the nodal- 
antinodal dichotomy observed in the pseudogap phase 
in underdoped cuprates. However, the anisotropy found 
here was still very small compared with the experimen- 
tal values observed in those materials. This is undoubt- 
edly related to our choice of initial interaction U = 3i. 
If we increase the on-site interaction of the model, we 
expect that this trend will be amplified further. Such a 
NFL character was observed for a finite doping range, af- 
ter which both nodal and antinodal quasiparticle weights 
became nonzero even when we extrapolated our RG flow 
to very low RG scales. This indicated that the result- 
ing FS, in this new phase, was fully reconstructed. With 
further doping, the system finally developed a d x 2_ y 2- 
wave pairing instability in the low-energy limit. This 
suggests the formation of (4,2 _ y 2 -wave superconducting 
ground state sufficiently away from half-filling. This lat- 
ter two-loop RG result therefore adds further support to 
the still much-debated existence of a d x 2_ y 2-wave super- 
conducting phase in the 2D Hubbard model and this is 
in qualitative agreement with the first results obtained 
in Refsi 12 ' 13 ! 14 within a one-loop RG scheme. 

Since our two-loop RG results have several analogies 
with the observed phenomena in the high-Tc supercon- 
ductors, a few remarks are in order. These materials are 
known to be Mott insulators at half-filling regime. This 
unambiguously implies that the strong-coupling regime 
should be an important ingredient to describe quantita- 
tively their observed properties. The main aim of this 
work however was not to achieve quantitative agreement 
with the experimental results. Our objective was to show 
that even the simplest version of the 2D Hubbard model 
is rich enough to capture some of the main features dis- 
played by those materials. In order to describe quantita- 
tively the experimental results of the high-Tc cuprates, 
it is essential to choose initial couplings larger than the 
ones considered in this work. Moreover, we should also 
take into account the effects produced by next nearest 
neighbor hoppings t' and t" . Nevertheless, it is reason- 
able to expect that, at least for the interaction strengths 
that apply to these compounds, the inclusion of t' and 
t" in the Hamiltonian should not alter qualitatively the 
low-energy dynamics of the system. In this sense, it is 
reasonable to suppose that the cuprates could be poten- 
tially associated with a universality class which is well 
captured by our model. As we have shown here, even our 
simple version of the 2D Hubbard model already displays 
some important and nontrivial aspects of the physics ex- 



hibited by these materials. 

On the other hand, there are still some open issues 
which were not addressed in this work. How the renor- 
malization of the FS induced by interactions affects the 
RG flow is still not considered in our RG scheme. In 
recent years, there has been some progress in analyz- 
ing this problem in quasi-one dimensional system s 33 ' 34 . 
However, a two-loop RG calculation of this effect in two- 
dimensional systems still remains a very difficult task 
and, to our knowledge, has not been calculated so far. 
Another important aspect which deserves future investi- 
gation concerns the actual importance of quantum fluc- 
tuations beyond the two-loop RG level for this model. 
Since three-loop RG calculations do not seem to be an 
easy prospect in the near future, we can nevetherless 
learn a lot about this question by comparing one-loop 
with two-loop RG calculations. One important result 
obtained here was the certainty that one-loop RG cal- 
culations in general tend to overestimate the effect of 
fluctuations near half-filling. As a result, our two-loop 
RG theory yields a better controlled approximation in 
the low-energy limit for this doping regime. In addition, 
our results suggest that the pseudogap phase might be 
indeed related to a strong-coupling regime in the model. 
However, like any other strong-coupling problem in cor- 
related systems, it is still not easy to have full analytical 
control in such case. In other words, there is still room 
for improvement before one can claim to have a complete 
theory for the strong-coupling regime of the 2D Hubbard 
model. However, given our encouraging results so far, 
we believe our two-loop RG approach could provide a 
good basis for a complementary view of this important 
problem starting from a weak-to-moderate coupling per- 
spective. 

Note added - Recently, we learned about a related 
study in the context of the 2D t — t' Hubbard model, 
RefM-. 
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